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The production of heavy quarkonia is a powerful tool to test our understanding of strong inter¬ 
action dynamics. It is well-known that the dominant production mechanism for heavy quarkonia 
with large transverse momentum is fragmentation. In this work we, analytically, calculate the QCD 
leading order contribution to the process-independent fragmentation functions (FFs) for a gluon to 
split into the vector (J/-i/>) and pseudoscalar (r] c ) S- wave charmonium states. The analyses of this 
paper differ in which we present, for the first time, an analytical form of the g —► J/if FF using a 
different approach (Suzuki’s model) in comparison with other results presented in literatures, where 
the Braaten’s scheme was used and the two-dimensional integrals were presented for the gluon FFs 
which must be evaluated numerically. The universal fragmentation probability for the g —> J/tf> is 
about 10~ 6 which is in good consistency with the result obtained in the Braaten’s model. 
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I. INTRODUCTION 

Heavy quarkonia, as the bound states of a heavy quark 
and antiquark, are the simplest particles when the strong 
interactions are concerned. Their production has a long 
history of theoretical calculations and experimental mea¬ 
surements [lj, especially, their production at high-energy 
colliders has been the subject of considerable interest dur¬ 
ing the past few years. Nevertheless, the production of 
heavy quarkonium is still puzzling us after almost forty 
years since the discovery of J/rp. New data have been 
taken at e + e _ , ep and pp colliders, and a wealth of fixed- 
target data also exist and with the advances in theory 
and tremendous amount of precise data from the Large 
Hadron Collider (LHC), it is an excellent time to study 
the physics of heavy quarkonium production. 

Nowadays, it is well-known that the dominant mecha¬ 
nism to produce the heavy quarkonia at high transverse 
momentum is fragmentation; the production of a parton 
with a large transverse momentum which subsequently 
decays to form a jet containing the expected hadron [2|. 
It is hence important to obtain the corresponding frag¬ 
mentation function (FF), in order to properly estimate 
the production rate of a specific quarkonium state. Be¬ 
cause of the simple internal structure of heavy quarkonia 
the perturbative QCD approximations to their FFs are 
well-defined in the nonrelativistic QCD (NRQCD) factor¬ 
ization framework 0]. To calculate the FFs, beside the 
current phenomenological approaches which are based on 
the x 2 analysis of experimental data (see our previous 
work [41]), there are two theoretical schemes which are 
based on the fact that the FFs for hadrons containing 
heavy quarks can be computed analytically using the 
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perturbative QCD (pQCD) S 0. In these theoreti¬ 
cal schemes the QCD improved parton model provides 
a great theoretical frame to extract the FFs. 

The first scheme for the inclusive production of heavy 
quarkonium has been developed by Bodwin, Braaten, 
and Lepage, who proposed a method on the basis of 
pQCD and the nonrelativistic quark model 0013. In 
this model, the FF of a heavy quark Q into the pseu¬ 
doscalar or vector heavy-light mesons Qq is defined as 
the cross section for producing a Qq -meson plus a light 
quark q with total four-momentum K^, divided by the 
cross section for producing an on-shell Q with the same 
three-momentum K, while K 0 —> oo. The fragmentation 
function in the Braaten’s model is defined as [7j 

i r vijw 2 

D{z,g 0 ) = J ds (1) 

where s = K 2 is the invariant mass of the meson, Mis the 
matrix element for producing mesons plus a light quark 
q , and M 0 is the matrix element for producing an on-shell 
Q. Here, D stands for the fragmentation function and z 
is the fragmentation parameter which refers to the lon¬ 
gitudinal momentum fraction of the quarkonium and go 
is a fragmentation scale. This model was applied to the 
fragmentation processes b —► B c and b B* in Ref. [ 8]. 
Another elaborate model is proposed by Suzuki 0, 110| 
which is based on the convenient Feynman diagrams and 
the wave function of the respective heavy meson, where 
the wave function includes the effect of long-distance in 
the fragmentation. In this model, the heavy FFs are cal¬ 
culated using a diagram similar to that in Fig. [3] so the 
analytical expression of FFs depends on the transverse 
momentum px of the parton which appears as a phe¬ 
nomenological parameter (e.g. see Eq. (12011 ). while in the 
Braaten’s model the integrations over all freedom degrees 
are performed. 

Note that, as soon as more than one hadron is appeared 
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in a hard scattering process, it is necessary to take into 
account the transverse momentum pt of the partons. For 
example, transverse momentum dependent FFs show up 
explicitly in several semi-inclusive cross sections, in par¬ 
ticular in azimuthal asymmetries. In the calculation of 
QCD corrections to these cross sections, the inclusion of 
the pt dependent FFs will be essential. Experimentally, 
based on the 1992-1993 run (run 1A), the CDF collabo¬ 
ration published data on their first measurement of the 
B-meson differential cross section da/dpx for the exclu¬ 
sive decays B + —> J/i/jK + and B° —> J/^K*° Ha. In 
jl2} . the pt distribution da/dpr is considered and the 
theoretical predictions are compared with the CDF data 
0 for which 5GeV< pr <20GeV. This is our main mo¬ 
tivation to study the transverse momentum dependent 
FFs in the Suzuki’s model. In Ref. using this model 
we computed an exact analytical expression of the frag¬ 
mentation function for c-quark to split into the D + /D° 
mesons to LO. There, we investigated that there is an 
excellent consistency between our result and the current 
well-known phenomenological models and also with the 
experimental data form BELLE and CLEO. 

In high energy processes, the main contribution of char- 
monium production results from gluon fragmentation, 
while the charm quark fragmentation contribution is 
much too small fl4 0- This point was confirmed by 
the comparison between the theoretical predictions and 
the CDF J/xjj production data. Therefore, in this work 
using the Suzuki’s model we focus on the gluon fragmen¬ 
tation into a vector ( J/ip ) and pseudoscalar ( iq c ) S-wave 
charmonium to leading order of perturbative QCD and 
we present, for the first time, an analytical expression 
for the transverse momentum dependent fragmentation 
functions of g —> r] c , J/^. 

Note that, in the past few years the g —»■ rj c , J/%j) FFs have 
been calculated numerically, using the Braaten’s model 
S0. In these papers, due to the lengthy and cumber¬ 
some expression of the g —> J/tf: FF the two-dimensional 
integrals have been presented that must be evaluated nu¬ 
merically. Moreover, since in the Braaten’s model inte¬ 
grations over all freedom degrees are performed, then the 
presented FFs are independent of the transverse momen¬ 
tum of the initial gluon. 

Our analytical expression for the g —»■ J/%1> FF will be 
compared with the numerical result presented in fl6| . 
Using the FF obtained, we also calculate the first two 
moments of FF which are of phenomenological interest 
and subject to experimental determination. They corre¬ 
spond to the g —> J/4’ branching fraction and the average 
energy fraction of the J/ip meson which receives from the 
gluon. 

This paper is organized as follows. In Sec. eh we explain 
our theoretical approach to calculate the FFs using the 
pQCD. Our analytical results of the g r] c ,J/ip FFs 
will be presented in the Suzuki’s model. In Sec. Mil our 
numerical results for the gluon FFs are presented. Our 
conclusion is summarized in Sec. EY1 



Figure 1: Feynman diagrams that contribute to charmonium 
production: (a) gg —> ccg at order-a^, (b) gg —> ccgg at 
order-Og Q- 


Figure 2: Feynman diagram of process qg —) qg* —¥ J/^gg 

Bi¬ 


ll. GLUON FRAGMENTATION INTO S-WAVE 
CHARMONIUM: PERTURBATIVE QCD 
SCHEME 

The theoretical schemes for calculating the charmo¬ 
nium FFs are based on the fact that the FFs for hadrons 
containing a heavy quark can be computed theoretically 
using perturbative QCD BSS0. Th ea-±r,c,J/il> 
FFs have been already calculated in Refs. , using 

the Braaten’s model where the fragmentation function is 
defined as in ©■ In Ref. Q, authors considered the typ¬ 
ical Feynman diagrams contributed to the production of 
charmonium states at the order (Fig. ©i) and at the 
order (Fig.©)). They have pointed out that in most 
regions of phase space, the virtual gluons in Fig.©) are off 
their mass shells by amounts of order pt , where pr refers 
to the large transverse momentum of charmonium state, 
and the contribution from this diagram is suppressed rel¬ 
ative to the diagram in Fig. ©i by a power of a s (pT) 
when the spin-singlet S-wave charmonium 7y c ( 1 /S'o) is con¬ 
sidered. Then, they have obtained the g g c FF at 
a s (fi = 2m c ) 2 -order. Authors have also calculated the 
FF for a gluon into J/if) considering the Feynman dia¬ 
gram for the process g* —>■ J/ipgg (Fig-© 1 ) to ct^. They 
have not presented an analytical form for the (z. p) 
and, instead, the result is given in a two-dimensional in¬ 
tegral form that must be evaluated numerically. In jl6j . 
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using the Braaten’s model authors have also obtained 
an integral form for the polarized and unpolarized initial 
FFs of gluon into the J/ip at a s (p = 2 m c ) 3 , considering 
a specific physical process shown in Fig. [5] 

Here, we focus on the gluon fragmentation into a point¬ 
like cc pair in the 1 So(rj c ) and 3 Si(J/ip) states and derive, 
for the first time, an analytical form of their transverse 
momentum dependent FFs, using the Suzuki’s model. 
The leading contribution to the short-distance scattering 
amplitude of the fragmentation process g —> J/ip arises 
from the partonic process g —> ccgg (Fig. [3]) and is of or¬ 
der a 3 , as its dominant decay process is J/ip —> 3 g. The 
hard scattering amplitude of the fragmentation g —>■ rj c 
results from the process g —> ccg (Fig. [4]) and is of order 

The fragmentation parameter z is normally defined in the 
Lorentz boost invariant form as z = ( E H +p |f)/ (E 9 +p 9 L ), 
which is known as the light-cone form. This definition 
is hard to be employed in the application of the gluon 
FF, because it involves the transverse momentum of the 
resulting heavy quarkonium. Instead, the following non¬ 
covariant definition is usually used approximately 


When the fragmenting gluon momentum |/c| —> oo, the 
definition © is equivalent to the light-cone definition 
0 - Therefore, we adopt the infinite momentum frame 
where the fragmentation parameter is defined as in 
In Ref. 0 ], authors analysed the uncertainties induced 
by different definitions of the momentum fraction z, in¬ 
cluding the covariant and non-covariant definitions and 
showed that the FFs corresponding to the light-cone def¬ 
inition of the fragmentation parameter £ are equivalent 
to the ones in the infinite momentum frame of gluon. 
Instead, the non-covariant definition © is used as an 
approximation unless |/c| —> oo. 

With the large heavy quark mass, the relative motion 
of the heavy quark pair inside the charmonium is effec¬ 
tively nonrelativistic [0 . For example, the squared rel¬ 
ative velocity of the heavy quark pair in the quarkonium 
rest frame is v 2 « 0.22 for the J/ip and v 2 ~ 0.1 for 
the T 0. The nonrelativistic assumption allows us to 
use a simple mesonic wave function which is the solution 
of the Schrodinger equation with a Coulomb potential 
|13| . A typical simple mesonic wave function is given 
in [20] which is the nonrelativistic limit of the Bethe- 
Salpeter equation with the QCD kernel. Here, according 
to the Lepage-Brodsky’s approach [2l[ we also neglect 
the relative motion of the heavy quark pair inside the 
charmonium states and we assume, for simplicity, that 
the quark pair are emitted collinearly with each other 
and move along the Z- axes. This assumption allows us 
to estimate the nonrelativistic mesonic wave function as 
a delta function form (more detail can be found in 0 )- 



Figure 3: Feynman diagram of process g —► J/ipgg. 


A. J/ip In gluon fragmentation 

To derive an analytical result of the function 
Dg^(z,g) which refers to the fragmentation of a gluon 
into the J/ip , we consider the Feynman diagram shown 
in Fig. [3l which is the subprocess of the physical pro¬ 
cess presented in Fig. [5] The spins (r^) and the four- 
momenta of meson and partons are also labelled in Fig. [3] 
According to our early assumption, the meson and its 
constituent quarks move along the Z-axes (fragmentation 
axes). Then we set the relevant four-momenta as 

s 'n = [ s 0i S T' s 'l\ t'g = [fo. *T) 4] 

Pn = [Po,Pt,Pl] Pg, = [P 0 ,0,P L }. (3) 

Considering the definition of fragmentation parameter, 
2 = E H /E 9 = Pq/po ©, we also may write the parton 
energies in terms of the initial gluon energy po as 

s 0 = Xizpo to = x 2 zp 0 
s'o = x 3 {l-z)po t' 0 = x 4 (l - z)p 0 , (4) 

where X\ = sq/Pq and x 2 = to/Po are the meson energy 
fractions carried by the constituent quarks which the con¬ 
dition x\+x 2 = 1 holds for them. Following Refs. |22l.l23|. 
in the nonrelativistic approximation we assume that the 
contribution of each constituent quark from the meson 
energy is proportional to its mass, i.e. X\ = m c /M and 
x 2 = iric/M where M = mj/^. Furthermore, we assume 
that the two gluon jets move almost in the same direction. 
This assumption is justified due to the fact that the very 
high momentum of the initial gluon is predominantly car¬ 
ried in the forward direction. Due to momentum conser¬ 
vation, the total transverse momentum of the two gluon 
jets will be identical to the transverse momentum of the 
initial gluon. Therefore, we have tip = = Pt/2. By 

this assumption one also has x 3 = X 4 = 1/2 in ©. 

In the Suzuki’s model the fragmentation function 
Dg^(z, p) at the initial scale po = rrij/^p, is obtained 
by squaring the total amplitude and integrating over fi- 
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nal state phase space fiol |23| 

S 

Xi 5 3 (P + s' +1' - p)d 3 Pd 3 s'd 3 t', (5) 

where, 77 is the spin of initial gluon and Tm is the prob¬ 
ability amplitude of the J/if) production. The amplitude 
Tm involves the hard scattering amplitude Th , which can 
be computed perturbatively from partonic subprocesses, 
and the process-independent distribution amplitude 
which contains the bound state nonperturbative dynamic 
of outgoing meson, i.e. 

T M (P,P,s',t') = f[dxi\T H (P,p, s',t',Xi)$M(xi,Q 2 ), 

( 6 ) 


Substituting Eqs. (0, 0 and (fllll) in 0 and carrying 
out the necessary integrations, the fragmentation func¬ 
tion reads 


Dg^iz,/! o) = -(ira s ) 3 (f M C F m%) 2 

f f EsiT^CP + s'+t'-p) 


(PoPos'ot'^iPo + s(, + t ’ 0 - Po ) 2 


d 3 Pd 3 s'd 3 t'. 

( 11 ) 


To obtain an analytical form of the fragmentation func¬ 
tion for gluon to split into the S-wave charmonia (i.e. 
rj c , J/i/j), we apply the scenario introduced in [26]. Ac¬ 
cording to this scenario if v(t) and u(s) are the Dirac 
spinors of the quarks forming the charmonium bound 
states, in the nonrelativistic approximation the projec¬ 
tion operator is defined as 


where, [dxi] = dx\dx 2 5{l — x\ — x 2 ) and Xi ’s are the 
momentum fractions carried by the constituent quarks. 
This scheme, introduced in |2(j l2ll|. is applied to absorb 
the soft behaviour of the bound state into the scatter¬ 
ing amplitude Th- The amplitude Th is, in essence, the 
partonic cross section to produce a heavy quark pair cc. 
with certain quantum numbers that, in the old fashioned 
perturbation theory is expressed as 

T = (47ra s (2w e ))imgC , F _ T _ 

2^/2P 0 p 0 s' 0 t' 0 {Po + t'o + Sq - Po) ’ 

where a s is the strong coupling constant, C F = v/5/12 is 
the color factor and T represents an appropriate combi¬ 
nation of the quark propagators and the spinorial parts 
of the amplitude. It reads 


r = GiG 2 


u{s, r 3 )f£>{</+ m c )f^v{t,r^) 


■ ( 8 ) 


Here, e/s are the gluon polarization vectors, G\ = 
1 /{q 2 — m 2 ) = 1 /( 2 s.s') and G 2 = l/{k 2 — m 2 ) = 1 /( 2 t.t') 
are proportional to the quark propagators. We put the 
dot products of the relevant four-vectors in the following 
form 


t.t' = s.s' 


■ 


4M(1 — z) 


Pt 


M (1 — z)m c 
4z 


(9) 


In 0, is the process-independent probability am¬ 
plitude to find quarks co-linear up to a scale Q 2 in the 
mesonic bound state, so that by working in the infinite- 
momentum frame it can be estimated as a delta function. 
Therefore, the distribution amplitude for a S-wave heavy 
meson with neglecting the Fermi motion, reads (2^, [25[ 


A s,s,(s,t) = v(t)u(s) oc (t/+m c )n s ,s z , (12) 

where nss 2 is the appropriate spin projection operator; 
noo = 75 for pseudoscalar state (r/ c ) and = ${S Z ) 

for vector state {J/ip)- The spin content of the polarized 
meson is then given by either 75 or (/, as could well be 
expected. This operator is convenient for our assump¬ 
tion in which we ignore the Fermi motion, so that the 
constituent quarks will fly together in parallel. There¬ 
fore, the spinorial part of the amplitude for formation of 
the vector charmonium state J/ij), which we denote by 
V, may be presented in the following form 

T v oc GiG 2 j((/+ rn c )^{(/ + m c )$i(l/ + m c )<fe j, (13) 

where q = s + s' and k = t + t' are the energy-momenta 
of the virtual intermediate quarks, e is the polarization 
four-vector of the meson J/ip which may be in a longi¬ 
tudinal state = e M (P, A = 0) or a transverse state 
e ( T )^ = e M (P, A = ±1). These components satisfy the 
relations; e(P, A).P = 0 , e^ T ^.P = 0 = x P and 

e(P, A).e*(P, A') = Therefore, for a vector char¬ 

monium with the four-momentum P M = [Po;0,P F ], the 
polarization four-vector is expressed as 

e {L) " = _L (j p i;0i0;j p o)) 

e (TV = -^(0;Tl,-i,0), (14) 

where Pl = sl + and Po = So + to- 
Now to obtain an analytical form of the J/%j) FF, in (fill) 
we perform a sum over the colors and the spins of gluons. 
Then the amplitude squared FF reads 


j. M * 2" i(x, - =j), do) 

where M is the meson mass and /m = i/12/M|\I/(0)| is 
the meson decay constant which is related to the non¬ 
relativistic mesonic S-wave function T(0) at the origin. 


^r y F y = 4GlGlTr[{1/+m c )(fy^{(/+m c )j''{I/+m c ) 

S 

x{l/+m c )'y v {q l +m c )j l j,<f{1/+m c )]. (15) 

Using the traditional trace technique, the trace may be 
expressed as the dot products of four-vectors. Here, we 










put the dot products of the relevant four-vectors in the 
following forms: 
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t .s = 0, 

p ' £t = ~71 ( } ’ 

M z 2 
= ~Tz + 2 M^' 

p.s = p.t = ^- + 


p.t' = p.s' = 


Pt 


f >/ r 1 / -i -\ 

= 'T-t 


• f J 

t.s = s.t = 


\ Z n Z o 

~2T mc + 80 ^) Pt ’ 


s'.e L = e L -t' = - 


M{ 1 - z) 


4 z 4M(l-z) 


Pt- ( 16 ) 


To proceed we need to specify the phase space integra¬ 
tions in ED- Note that 


/ 


d 3 P5 3 (P+t'+ s'^p) 
Po(Po + s o + — Po) 2 


Po 


(M 2 + 2p.t' + 2 p.s' - 2 s’.t’) 2 ‘ 


(17) 


Here, instead of performing the transverse momentum 
integrations we replace the integration variable by its av¬ 
erage value (pt) in each case, which is a free parameter 
and can be specified experimentally. Therefore we can 
write 

J F{z,s' T )d 3 s' = j F{z, s' T )ds' L d 2 s' T 

« m 2 c s' 0 F(z,(s'S)) = m 2 c s' 0 F(z,^ {p 2 T )), 

(18) 

and 

J F'{z,t' T )d 3 t' = J F'{z,t' T )dt' L d 2 t' T 

~ m 2 t' 0 F'(z, (t'S)) = m%F'{z , ^ (p 2 T )). 

(19) 

Finally, putting all in (fill) and by assuming M = 2 m c in 
the nonrelativistic limit, we obtain the longitudinal and 
the transverse components of the g —>■ J/ip FF, as 


Pg^tj!^ z i A*o) 


8 Niza^ 

9(z,(Pt)) 


z( 1-z) 


- 16- 


.1 -z 2 


(2z 3 + 4z 2 + Zz - 4)m A c p 2 T 


{1-z) 3 


z{17z 2 — 25 z + 8)m 2 c p T — z 4 p? 


( 20 ) 



Figure 4: Feynman diagram of process g —t g c at af-order. 


and 


P g—P®) 


9(z,(Pt)) 

16 


- 32- 


1 — z 2 


{1-z) 3 


+- -(3 z 2 -\- 2 z 2 )ui*/pj. T 

1 — z 


2z 2 {—4z 2 — 7 z + 11 )m 2 c PT + 3 z A p^ 


( 21 ) 


where N\ = (7r 3 m®/| f C'J.)/216, and 
9(z,(pt)) = 


zm c 2 _j_ M{1 — z)m c 


2M{1 — z) 


M 2 


2 z 


1 2 


1 - 2 


Pt 


( 22 ) 


Note that, the fragmentation function for a vector char- 
monium J /ip is the sum of the longitudinal and twice the 
transverse components, i.e. Dg^{z,p 0 ) = 2+ Djp. 
Finally, we point out that to impose the effect of quarko- 
nium spin into the calculation, a second scenario is de¬ 
fined in ( 22 I . [27l | , where in the nonrelativistic approxima¬ 
tion the spin projection operator is expressed as 

A 5A (P) = i|(f + M)H SSz , (23) 

where, /m is the meson decay constant and M and P are 
the mass and the four-momentum of meson bound state, 
respectively. In the nonrelativistic limit {p = 21/M = 
2 m c ), these two scenarios El ED are the same and, in 
conclusion, the sum of transverse and longitudinal polar¬ 
isations of the J/ip should be identical. 


B. r/ c In gluon fragmentation 

Following Ref. Q, the fragmentation function for a 
gluon to split into the pseudoscalar S-wave charmonium 
r] c can be calculated by considering Fig. [l]i. To derive 
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an analytical result of the g —> r) c FF we consider the 
Feynman diagram shown in Fig. [H where a gluon forms 
a bound state cc with a gluon produced through a sin¬ 
gle c-quark. The four-momenta of meson and partons 
are also labelled. According to our previous assumption 
the meson and its constituent quarks move along the Z- 
axes (fragmentation axes). Then one can set the relevant 
four-moment a as 

Pn = \Po, Pt ,Pl\ kp = [k Q , p T ,k L ] 

PP = bo. 0 ,Pl] K = k 'h\ 

P lt = [P 0 ,O,P L ]. (24) 

Considering the definition of fragmentation parameter 
©, we may write the parton energies in terms of the 
gluon energy p 0 as; p ' 0 = x\zp 0l k ' 0 = x 2 zp 0 , k 0 = (l-z)p 0 
where X\ and x 2 are the meson energy fractions carried 
by the constituent quarks. Following our early discus¬ 
sion, the contribution of each constituent quark from the 
meson energy is proportional to its mass, i.e. X\ = m c /M 
and x 2 = ms/M where M = mj/^. 

We start with the definition of the fragmentation func¬ 
tion as 

= TT2T^ / l T M| 2 <5 3 (P + k-p)dWk, (25) 

Q s 

where the probability amplitude Tm is related to the hard 
scattering amplitude Th and the distribution amplitude 
as in ©. In view of our early discussion in this 
section, we propose a delta function as m for the am¬ 
plitude 4 >m- 

Using the perturbation theory, the amplitude Th is writ¬ 
ten in the following form 

T ^ 4tt a s {2m c )m 2 c C F _ T 

2\/2P 0 p 0 k 0 (Po + ko - Po) 

where T represents the spinorial parts of the amplitude 
as 


F = Gs u(fc , ,r 2 )^(fc)(g / +m c )^i(p)n(p',ri) b (27) 


part of the amplitude (1571) is presented in the following 
form 


r p 



{¥ 


m c h 5 <&(k)(</ + m c )q{(p) >• 


(29) 


By performing a sum over the colors and the spins of 
gluons, the amplitude squared reads 


^r p r p = 128G 2 m 2 d 2m 2 + 2 k.k'+p'.(k + k') [.(30) 


Next we consider the phase space integrations (1281) . Fol¬ 
lowing our previous approach, we have 


d 3 P< 5 3 (P + k — p) 
Po(Po + k 0 -po ) 2 


Po 

(M 2 + 2 p.k ) 21 


(31) 


and 


J F(z,p F )d 3 k = m 2 c k 0 F(z, (p%)). (32) 


In order to get the correct result for the fragmentation 
function of the process <7 —>• 77 c it is necessary to consider 
a second diagram that can be obtained from Fig. [4] by 
interchanging the two vertices where the gluons attached 
to the heavy quark lines. 

In conclusion, in the referred scenario the fragmentation 
function of the process <7 —>■ 77 ,- is expressed as 


A?-»?c( 2 >Mo) = (192 to c ) 


30:(Pt» 


2 m r 


z Pt 


M (1 — z) 


M( 1 - z) 


(33) 


where N 2 = (irm 3 f^Cp ) 2 /864, and 


9(z,(Pt )) = 


(M 2 


z 2 Pt z 2 Pt + AF 2 (1 - zf 


1 - 2 


)( 


Mz(l-z) 


(34) 


Here, G = l/(q 2 — to 2 ) = 1/(2 k.k') is proportional to the 
quark propagator and C F = 1 /(2v^3) is the color factor. 
Putting all in ((55l) . the fragmentation function reads 

Dg c (z,P 0 ) = 1 ( 7 Ta s f M C F m 2 c ) 2 

x /'i^££rf 3 k f ^(P + F-P) c/3 p 
J Poko J Po(Po + k 0 - Po r 

(28) 

Considering the scenario introduced in previous section 
urn where the spin projection operator for a pseu¬ 
doscalar charmonium state r] c is Hqo = 75, the spinorial 


III. RESULTS AND DISCUSSION 

We are now in a position to present our phenomenolog¬ 
ical predictions for the gluon fragmentation into the rj c 
and J /ip, by performing a numerical analysis. In general, 
the fragmentation function (z, p) depends on the fac¬ 
torization scale p. Here, we have set the scale in the FF 
and in the running coupling constant ct s (p) to p = 2 m c 
which is the minimum value of the invariant mass of the 
fragmenting gluon. Therefore, the functions (l 20 l [ 2 TT) and 
(1551) should be regarded as models for the g g c ,J/ip 
FFs at the initial scale p 0 = 2m c . For values of p much 
larger than po, the initial FFs should be evolved from the 
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Figure 5: The fragmentation function for the process g —> 
r/c, considering two values of the transverse momentum pr = 
10,15 GeV. The initial scale is set as yo = 2 m c . 


Figure 6 : The fragmentation function for the process g —> J/ip 
at 0 's. Here, D 3 J^ (z,y o) = 2Z)J + Dg where the transverse 
component Z)J is given in (1201) and the longitudinal one Dg 
is given in (HD- 


scale yo to the higher scale y using the Altarelli-Parisi 
equations |28l. [ 2 !)i| 


d 

-+h(z,h) = Y, 


dy 


Pi^j{z/y,y)Dj^ H (y,n), 


(35) 


where Pi^j(x, y) are the Altarelli-Parisi functions for the 
splitting of the parton of type i into a parton of type j 
with momentum fraction x. The only boundary condi¬ 
tion on this evolution equation is the fragmentation func¬ 
tion D^h {z, yo) at the scale yo = 2 m c . 

For numerical results, we take m c = 1.5 GeV, mj= 
3096.9 MeV, m Vc = 2983.6 MeV, a s (2m c ) = 0.26 and 
/m(cc) = 0.48 GeV 0], 

In Fig. [5j our prediction for the <7 —M j c FF (1551) by con¬ 
sidering two values of the gluon transverse momentum 
(pr = 10,15 GeV) is shown. It shows that the frag¬ 
mentation function distribution relies on the momentum 
Pt of the initial gluon. In Fig. [ 6 l the behaviour of the 
g —> J/ij) FF is studied. Here, Dg^ is the convenient 
summation of the longitudinal and transverse fragmenta¬ 
tion functions as D J g ^{z, yo) = 2 Dj + Dg, see Eqs. (l20l 
EH)- Our result shown in Fig. E3 is in acceptable agree¬ 
ment with the result presented in Fig. 3 of Ref. [16f, when 
the non-covariant definition of fragmentation parameter 
© is applied. In both results, the peak position of the 
FF occurs at 0.22 when pp = 10 GeV is considered, 
and the maximum value of the FF is D'!/' h m 0.8 x 10 -5 . 
Besides the g —> J/-0 FF itself, also its first two moments 
are of phenomenological interest. They correspond to the 
g —> J/%1) branching fraction 

B(y) = [ dzD g ^j/^{z,y), (36) 

Jo 


and the average energy fraction that the J/tp meson re¬ 
ceives from the gluon 

1 r 1 

( z ) ( m ) = zdzD g ^j^{z,y). ( 37 ) 

Indeed, an order of magnitude estimate of the gluon 
fragmentation contribution to J/tj) production in any 
high transverse momentum process can be obtained by 
multiplying the cross section for producing gluons with 
Pt > 2 m c by the branching fraction, which refers to the 
fragmentation probability. Our result for the g —> J/ip 
branching fraction is Bj/^,(2m c ) = 2.94 x 10 -6 and 
(z) j/tJj (2 m c ) = 0.277 which this initial fragmentation 
probability can be compared with the result presented 
in [2J where B(2m c ) = 3.2 x 10 -6 . 

Our results may be directly applied to the S-wave 
bottomonium sates T( 3 Si) and 77 b( 1 5 * 0 ), except that m c 
is replaced by mb = 4.5 GeV and the decay constant 
/m( 66 ) = 0.33 GeV (3C)1 | is the appropriate constant for 
the bottomonium states. Since the b-quark is heavier 
than the c-quark, we expect that the peaks of the frag¬ 
mentation functions shift significantly toward higher val¬ 
ues of z. Our results for the g —> 77 c ,«//■*/» FFs will be 
checked by the comparison between the theoretical pre¬ 
dictions and the experimental measurements of the heavy 
meson cross sections at the LHC. 


IV. CONCLUSION 

Understanding hadronization, the process by which a 
parton evolves into a hadron, is complicated by its in¬ 
trinsically nonperturbative nature. In hadron colliders, 
at sufficiently large transverse momentum of the heavy 
quarkonium production the direct production schemes 










are normally suppressed while the fragmentation mech¬ 
anism becomes dominant. The fragmentation refers to 
the process of a parton with high transverse momentum 
which subsequently decays into the expected hadron Q. 
Beside the phenomenological approaches, it is well-known 
that the fragmentation function which describes this pro¬ 
cess can be calculated using perturbative QCD. In this 
paper we, for the first time, gave out an analytical form 
for the leading color-singlet contribution to the fragmen¬ 
tation function for a gluon to split into the vector and 
pseudoscaler S-wave charmonia at the initial 

scale fj, = 2 m c . We used a different model in getting 
them from the Braaten’s model applied in other litera¬ 
tures, see Refs. @,13. Our results depend on the trans¬ 
verse momentum of the initial gluon and shows a sim¬ 


ple analytical form whereas in the Braaten’s model, the 
integrations over all freedom degrees are performed and 
due to the lengthy and cumbersome expression of the FFs 
the results are presented as the two-dimensional integrals 
that must be evaluated numerically. Since the transverse 
momentum dependent FFs show up explicitly in semi- 
inclusive cross sections, therefore in the QCD corrections 
the inclusion of these dependent FFs will be necessary. 
Our result for the g —► J/ip FF is in acceptable consis¬ 
tency with the numerical result presented in Ref. 3, 
when one uses the normal definition of the fragmenta¬ 
tion parameter We also found that the fragmenta¬ 
tion probability of a high energy gluon splitting into J/tj) 
is about 10~ 6 . 
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